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ABSTRACT 



We present a calculation of the rare decay mode — > 77, in which the long distance contri- 
■ butions are expected to be dominant. Using the Heavy Quark Chiral Lagrangian with a strong g 

^] , coupling as recently determined by CLEO from the D* — > Dtt width, we consider both the anomaly 

I contribution which relates to the annihilation part of the weak Lagrangian and the one- loop vr, K 

diagrams. The loop contributions which are proportional to g and contain the oi Wilson coefficient 
are found to dominate the decay amplitude, which turns out to be mainly parity violating. The 
1-^ I branching ratio is then calculated to be (1.0 ± 0.5) x 10~^. Observation of an order of magnitude 

larger branching ratio could be indicative of new physics. 



PACS number(s): 13.25.Ft, 14.70.Bh, 12.39.Hg, 12.39.Fe 



1 Introduction 



With the new data coming and expected from the B-factories, there is very strong emphasis and 
activity in the field of B-physics in all its aspects. This includes the rare decays of B-mesons, which 
are considered as an attractive source for possible signals of new physics. In contradistinction to 
the growing efforts to understand the decay mechanisms of rare B decays, studies of rare D decays 
have received less attention in recent years. Partially this is because theoretical investigations of D 
weak decays are rather difficult, also due to the presence of many resonances close to this energy 
region. The penguin effects on the other hand, which are very important in B and also in K decays, 
are usually suppressed in the case of charm mesons due to the presence of d, s, b quarks in the loop 
with the respective values of CKM elements. 

Nevertheless, D meson physics has produced some interesting results in the past year. Experi- 
mental results on time dependent decay rates of K^-k' by CLEO and — > K~^K~ and 
K~7r~^ by FOCUS have stimulated several studies on the Z?*^ — Z)*^ oscillations |^. The 
recently measured D* decay width by CLEO has provided the long expected information on the 
value of D*Dir coupling. Among the rare D decays, the decays D Vj and D V{P)l^l~ are 
subjects of CLEO and FERMILAB searches Q. On the theoretical side, these rare decays of charm 
mesons into light vector meson and photon or lepton pair have been considered lately by several 
authors (see, e.g., p-ljlll and for radiative leptonic D meson decay see |jl^). The investigations 
of D — > Vj showed that certain branching ratios can be as large as 10~^, like for 
Df /9+7 [P, However, the decays which are of some relevance to the — > 27 mode studied 
here, like — > p^'-f, — > loj, are expected with branching ratios in the 10~^ range |13]. Thus, it 
is hard to believe that the branching ratio of the 27 decay mode can be as high as 10~^ in 
the Standard Model (SM), as found by [14|. Apart from this estimation, there is no other detailed 
work on — > 27 in the literature, to the best of our knowledge. 

On the other hand, in the B and K meson systems there are numerous studies of their decays 
to two photons. For example, the Bs — > 77 decay has been studied with various approaches within 
SM and beyond. In SM, the short distance (SD) contribution leads to a branching ratio 
B{Bs — > 77) — 3.8 X 10^''. The QCD corrections enhance this rate to 5 x 10^^ [16]. On the other 



hand, in some of the SM extensions the branching ratio can be considerably larger. The two Higgs 
doublet scenario, for example, could enhance this branching ratio by an order of magnitude |17|. 
Such "new physics" effects could at least in principle be dwarfed by long distance (LD) effects. 
However, existing calculations show that these are not larger than the SD contribution [|l8i, which 



is typical of the situation in radiative B decays |19|. In the system the situation is rather 
different. Here, the SD contribution is too small to account for the observed rates of Kg — > 27, 
Ki — > 27 by factors of ~ 3 — 5 [20 1, although it could be of relevance in the mechanism of CP- 
violation. Many detailed calculations of these processes have been performed over the years (see 
recent refs. [^0|-|2^ and refs. therein), especially using the chiral approach to account for the pole 
diagrams and the loops. These LD contributions lead to rates which are compatible with existing 
measurements. 

Motivated by the experimental efforts to observe rare D meson decays, as well as by the lack of 
detailed theoretical treatments, we undertook an investigation of the 77 decay. The short 

distance contribution is expected to be rather small, as already encountered in the one photon 
decays |6|, 0, hence the main contribution would come from long distance interactions. In order 
to treat the long distance contributions, we use the heavy quark effective theory combined with 
chiral perturbation theory (HQxPT) ||2^. This approach was used before for treating D* strong 
and electromagnetic decays |p5|-[p^]. The leptonic and semileptonic decays of D meson were also 
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treated within the same framework(see [^] and references therein). 

The approach of HQxPT introduces several coupling constants that have to be determined from 
experiment. The recent measurement of the D* decay width Q has determined the D*Dtt coupling, 
which is related to g, the basic strong coupling of the Lagrangian. There is more ambiguity, however, 
concerning the value of the anomalous electromagnetic coupling, which is responsible for the D*D^ 
decays 27], as we shall discuss later. 



Let us address now some issues concerning the theoretical framework used in our treatment. For 
the weak vertex we shall use the factorization of weak currents with nonfactorizable contributions 
coming from chiral loops. The typical energy of intermediate pseudoscalar mesons is of order m£)/2, 
so that the chiral expansion p/A^ (for > 1 GeV) is rather close to unity. Thus, for the decay 
under study here we extend the possible range of applicability of the chiral expansion of HQxPT, 
compared to previous treatments like D* — > Dtt, D* — > D-f [^] or D* — > D'yy |27|, in which a 
heavy meson appears in the final state, making the use of chiral perturbation theory rather natural. 
The suitability of our undertaking here must be confronted with experiment, and possibly other 
theoretical approaches. 

At this point we also remark that the contribution of the order 0{p) does not exist in the 
£)0 _^ decay, and the amplitude starts with contribution of the order 0{p^). At this order the 
amplitude receives an annihilation type contribution proportional to the 02 Wilson coefficient, with 
the Wess-Zumino anomalous term coupling light pseudoscalars to two photons. As we will show, 
the total amplitude is dominated by terms proportional to ai that contribute only through loops 
with Goldstone bosons. Loop contributions proportional to 02 vanish at this order. We point out 
that any other model which does not involve intermediate charged states cannot give this kind of 
contribution. Therefore, the chiral loops naturally include effects of intermediate meson exchange. 
The chiral loops of order 0{p^) are finite, as they are in the similar case of iiT — > 77 decays 
23]. The next to leading terms might be almost of the same order of magnitude compared 



to the leading 0{p^) term, the expected suppression being approximately p^/A^. The inclusion 
of next order terms in the chiral expansion is not straightforward in the present approach. We 
include, however, terms which contain the anomalous electromagnetic coupling, and appear as next 
to leading order terms in the chiral expansion, in view of their potentially large contribution (as 
in B*{D*) —>■ B{D)jj decays considered in ]|2^). As it turns out, these terms are suppressed 
compared to the leading loop effects, which at least partially justifies the use of HQ^PT for the 
decay under consideration. Contributions of the same order could arise from light resonances like 
p, K* , ao(980), /o(975). Such resonances are sometimes treated with hidden gauge symmetry (see, 
e.g., [^]), which is not compatible with chiral perturbation symmetry. Therefore, a consistent 
calculation of these terms is beyond our scheme and we disregard their possible effect. 

Our paper is organized as follows: in Section ^ we present the basic features of the model. We 
give the results and their discussion in Section ^ and conclude with a summary in Section 



2 The theoretical framework 



The invariant amplitude for D'^ — > 77 decay can be written using gauge and Lorentz invariance in 
the following form: 



M 



^2 ^1 

ki-k2' 



eifj.e2u, 



(1) 



where ^ is a parity violating and M^^^ a parity conserving part of the amplitude, while ki(^2)^ 
ei(2) are respectively the four momenta and the polarization vectors of the outgoing photons. 
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In the discussion of weak radiative decays 4 — > 577 or g' — > 57 decays, usually the short (SD) 
and long distance (LD) contribution are separated. The SD contribution in these transitions is 
a result of the penguin-like transition, while the long distance contribution arises in particular 
pseudoscalar meson decay as a result of the nonleptonic four quark weak Lagrangian, when the 
photon is emitted from the quark legs. Here we follow this classification. In the case of 6 ^ 577 
decay it was noticed that without QCD corrections the rate V{b — > s77)/r(6 — > 57) is about 
10~^. One expect that a similar effect will show up in the case of c — > U77 decays. Namely, 



according to the result of |28| the largest contribution to c — > M77 amplitude would arise from the 
photon emitted either from c or u quark legs in the case of the penguin-like transition c — > U7. 
Without QCD corrections the branching ratio for c ^ n7 is rather suppressed, being of the order 
10"^'' The QCD corrections [^] enhance it up to order of 10~^. 

In our approach we include the c — > U7 short distance contribution by using the Lagrangian 

L = -^VusV:,C%f ^^F,.m. {ua>^''\{l + 75)0) , (2) 

where rric is a charm quark mass. In our analysis we follow ^] and we take Cj!^f = (-0.7 + 
2i) X 10-2. 

The main LD contribution will arise from the effective four quark nonleptonic AC = 1 weak 
Lagrangian given by 

^ = "Tf E VuqV:^[ai{qT^'c){uT^q) + a2(nr^c)(gr^(z)] , (3) 

q=d,s 

where T'^ = 7^(1 — 75), aj are effective Wilson coefficients [^], and Vg-q. are CKM matrix elements. 
At this point it is worth pointing out that long distance interactions will contribute only if the 
SU{3) flavor symmetry is broken, i.e. if 7^ niii. Namely, due to VudVcd ~ — KjsK*s) if ^rf = "^s 
the contributions arising from the weak Lagrangian (^) cancel. 

Now, we turn to describe some of the basic features of the HQxPT. This model would serve 
us as a hadronized counterpart of the quark effective weak Lagrangian. One has the usual 0{p^) 
chiral Lagrangian for the light pseudoscalar mesons 



(4) 



where S = exp(2in//) with 11 = containing the Goldstone bosons 7r,K,ri, while the 

trace tr runs over flavor indices and. — d.ia;g(?77,^, md,ms) is the current quark mass matrix. 
From this Lagrangian, we can deduce the light weak current of the order 0{p) 

= -i^tr(Sa^StA^), (5) 

corresponding to the quark current = qLlfj.^^ QL- (A^ is an SU(3) flavor matrix.) 

In the heavy meson sector interacting with light mesons we have the following lowest order 0{p) 
chiral Lagrangian 

££ = - Tv{Haiv-DabHb) + g TT{HaHt-f^Al 75) , (6) 

where D'^f^Hf, = d^Ha — HhVl^^, while the trace Tr runs over Dirac indices. Note that in (^) and the 
rest of this section a and b are flavor indices. 
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The vector and axial vector fields and in (|6|) are given by: 

V, = liCd,^^ + c^d.o A, = ^(e^a.e - m^) > (7) 

where ^ = exp (ill//). The heavy meson field Ha contains a spin zero and spin one boson 

Ha =P+{P^,al^ - P5al5), (8) 
Ha =7°(^a)^7° = [Pla^ + ^aTs] P+ , (9) 

with P± = (1 lb 7^?j^)/2 being the projection operators. The field ^5(^5) annihilates (creates) a 
pseudoscalar meson with a heavy quark having velocity v, and similar for spin one mesons. 
For a decaying heavy quark, the weak current is given by 

Ja = Qal^LQ , (10) 

where L = (1 — 75)72 and Q is the heavy quark field in the full theory, in our case a c-quark, and 
q is the light quark field. 



From symmetry grounds, the heavy-light weak current is bosonized in the following way |24| 



J^ = '-^Tr[^^LH,cl], (11) 
where a is related to the physical decay constant fo through the well known matrix element 

(0|n7S5c|Z)°) = -2(0| J,^|L>°) = imDV^fn • (12) 



Note that the current ( [111) is 0{p^) in the chiral counting. 

In the calculation of short distance contribution (^) there appears the operator (^tc^i/| (1+75)0). 
Using heavy quark symmetry this operator can be translated into an operator containing meson 
fields only [36| 



|(l + 75)c) - ^Tr[a^4(l+75)/f4j. (13) 



[ua 
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The photon couplings are obtained by gauging the Lagrangians (^ , (^ and the light current (|5|) 
with the U{1) photon field B^. The covariant derivatives are then D'^f^Hb = d^Ha + ieB^{Q'H — 
HQ)a - HbVf^^ and = + ieBfj_[Q,^] with Q = diag(|, -i, -i) and Q' = | (for our case 
of c quark). The vector and axial vector fields (^) change after gauging and now they read = 
^(^^AtC^ + C^^mO ^J^d = §(^^^^/i^ — ^T^/iC^)- The light weak current (|5|) contains after gauging 
the covariant derivative instead of d^. However, the gauging procedure alone does not introduce 
a coupling between heavy vector and pseudoscalar meson fields and the photon without emission or 
absorption of additional Goldstone boson, which is needed to account, for example, for D* — > D'j. 



To describe this electromagnetic interaction we follow [26| introducing an additional gauge invariant 



contact term with an unknown coupling /3 of dimension -1. 

£^ = T¥ HaHba'^-'F^^Ql - ^Q' Tr Haa^^'HaF^,, (14) 

where = ^iC^QC + CQ^^) and F^^ = d^B^ - d^B^. The first term concerns the contribution of 
the light quarks in the heavy meson and the second term describes emission of a photon from the 
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Figure 1: One loop diagrams (not containing /3-like terms (|l^)) that give vanishing contributions. 
The dashed hne represents charged Goldstone bosons flowing in the loop {K^ , 7r~^), while the double 
line represents heavy mesons, D and D* . 



heavy quark. Its coefficient is fixed by heavy quark symmetry. From this "anomalous" interaction, 
both H*Hj and H*H*'^ interaction terms arise. Even though the Lagrangian ([T^ ) is formally 
l/mq ~ ruq suppressed, we do not neglect it, as it has been found that it gives a sizable contribution 
to D*{B*) D{B)jj decays [^]. In the case of — > 77 it gives largest contribution to the 
parity conserving part of the amplitude, however it does not contribute to the decay rate by more 
than 10%, as will be shown later. The Lagrangian (jlj) in principle receives a number of other 
contributions at the order of l/mq, however, these can be absorbed in the definition of /? for the 
processes considered [^]. 



3 Results 



The decay width for the 



77 decay can be obtained using the amplitude decomposition in 



^(|M<-)|^ + i|M<-.^„i), 



(15) 



The short distance contribution to 77 decay width is estimated using c — > M7 effective 

Lagrangian (|2|), ( p!3| ) with one photon emitted from leg via term (0). The parity violating 
part of the short distance amplitude is 



3/2 ^ 

while the parity conserving part of the amplitude is 



l + 2A*/mL>- 



127r2 ^ 



niD + 2A* 



(16) 



(17) 



where A* = m^jo* —rrij^o. Turning now to the long distance contributions, we depict in Figs. |^ and 
§ the loop diagrams arising to leading order 0{p^) by using Eqs.(§)-(P) and (11). The circled cross 
indicates the weak interaction. In Figure || we grouped all diagrams which vanish by symmetry con- 
siderations. All nonvanishing contributions are assembled in Fig. |2[ We denote the gauge invariant 
sums corresponding to nonvanishing diagrams of Fig. || by M-^^ = M-^^ (the gauge invariant 
sums are sums of diagrams in each row of Fig. 0), where +(— ) denotes parity conserving (violating) 
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Figure 2: One loop diagrams, not containing beta-like terms (14), that give nonvanishing contri- 
butions to the 77 decay amplitude. Each sum of the amplitudes on diagrams in one row 
Mi = ^ • Mij is gauge invariant and finite. Numerical values are listed in Table |^. 




Figure 3: Anomalous contributions to 77 decay. The intermediate pseudoscalar mesons 

propagating from the weak vertex are TT^,r],r]'. 
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part of the amplitude, as in (|T5|). The parity violating sums, which arise from the first term in @ are 



; Gf 2 1 



— A7^2(mi,, ^ + a:) - 2G3(m^, + a:, -^^ 

-Gf o 1 r 



(18) 



)) + (19) 



^7-7r7^2(m^, ^ + AS) - 2G3(m^, mo + A^, -H^)) 

+ An ^ 



G 1 



(20) 



with 



777^ r 1 
/(m,A, ms) = Go(m, A + ^, HLa ) _ Go(m, A, ^ 



Go(m, A + I^, HLa) _ AGo(m, A, I^)] + ^ + ^ 

— A^o("i,m|)) + (^ + — ) Go(m,A + ^,mi5) +m2Mo("i,A + ^,mD) 



+ 



(m2-A2)r 1 



+ (a - ^)M,{m, -\{^- A)iVo(m,ml,) 

(2A + mz)) ^ , (3m?./2 + 3AmD + 2A2-2m2) . ^ 

4m£)(A + m£)) 2m|j(mD + 2A) ^ 

(mi^A - 2m2 + 2A2) 
H ■ 2 A ^2("^, A). 



The parity conserving parts of the amplitude M-~^^ vanish for the diagrams on Fig. y. We denoted 
Ag*^ = m (*) —m^o, while functions l2{m, A), GQ{m, A, v-fc), Gs{m, A, f-A;), N(){m, k"^), Go(m, A, w- 

A;), Mo(m, A, A;), M2{m,ki-k2), M4(m, fci- ^2) are presented in the Appendix^. Note that the 
sums of amplitudes (p!^)-(pO|) are gauge invariant and finite. This is expected, since one cannot 
generate counterterms at this order. There is no ^ dependence apart from the one hidden in ai 
even though jj, appears in the above functions, but it cancels out completely. Note also, that the 
one loop chiral corrections vanish in the exact SU{3) limit, i.e. when rriK — > m^, as it is expected. 
One should note that taking the chiral limit (i.e. ms^md 0) is not unambiguous. Namely, in the 
combined heavy quark effective theory and the chiral perturbation theory, beside chiral logarithms 
there are also functions of the form F{mq/A) whose value depend on the way one takes the limit 
(see e.g. Ref. ||). 

We remark that there exist additional diagrams of the same order in the chiral expansion as 
ones given on Fig. ^, but proportional to 02 part of the effective weak Lagrangian (^). In these 
additional diagrams the chiral loop is attached to the light current in the factorized vertex, while 
the photons are emitted from the pseudoscalars in the loop, or they come from the weak vertex. 
However, the amplitudes of these diagrams vanish due to Lorentz symmetry. 



(21) 
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The contribution coming from the anomalous coupUng tt^jj, 7777, ??'77 (Fig. Isl) is 



M 



(+) 



Anom. 



, Gf e2 



E 



rriD 



P=iT'>,ri,rj' ^ 



nil 



-Kp 



(22) 



COS © 'I 



* / sine I _\/2cose\l r_x/2cose _ 4sine i 



+ VusV. 



T^l r \r T/*/'sin0 I cos©"! I T A* / sin © cos © N 1 f \/2 sin © i 



V3 



V3 



V3 



where 9 = —20° it 5° is the r] — rj' mixing angle and we have set fj^ 



parameters reproduces the experimental results for the vr^ 



77, r] 



fv8 ~ 

^ 77, 



' fvo- 
and T]' 



4 cos e 1 

This choice of 



77 decay 



width [34|. In the numerical evaluation we use the values of a and g obtained within the same 

The coupling g is extracted from existing experimental data 
Recently CLEO Collaboration has obtained the first measurement of D*~^ decay 



framework as in [25, 
on D* Dtt. 

width r(L>*+) = 96 ± 4 ± 22 keV [|] by studying the D* L>°7r+. Using the value of decay 
width together with branching ratio Br{D*~^ D^ir^) = (67.7 it 0.5)% one immediately finds at 
tree level that g = 0.59 it 0.08. The chiral corrections to this coupling were found to contribute 
about 10% |25, In order to obtain the a coupling, we use present experimental data on Ds 
leptonic decays. Namely, at the tree level there is a relation fu = f^^ = a/y/niD- From the 
experimental branching ratio Dg — > //f^x and the Dg decay width [^4| one gets /d^ = 0.23 it 0.05 
GeV and a = 0.31 it 0.04 GeV'^''^. The SU{3) breaking effects in the form of chiral loops and 
the counterterms can change the extracted value of a. One chiral loop corrections can amount 
to about 40% when g is taken to be 0.59. This value might be changed by the finite part of the 
counterterms. However, the contributions coming from counterterms are not known and due to the 
lack of experimental data they cannot be fixed yet. In our calculation we take a = 0.31 GeV3/^ 
keeping in mind that the chiral corrections might be important (for instance setting counterterms 
to zero in one loop calculation one gets a = 0.21 ± 0.04 GeV^/^ using g = 0.59, for details see 

. We present 



Appendix B of |32|). For the Wilson coefficients ai we take 1.26 and a2 
the numerical results for the one loop amplitudes in Table ||. 



-0.47 





Ml ^ [xlO-i° GeV] 


m|+^ [x 10-10 GeV-i] 


Anom. 





-0.53 


SD 


-0.27 -0.81f 


-0.16 -0.47i 


1 


3.55 +9.36i 





2 


1.67 





3 


-0.54 +2.84i 







4.41 +11.39i 


-0.69 -0.47i 



Table 1: Table of the nonvanishing finite ampfitudes. The ampfitudes coming from the anomalous and short distance 
{Cjl^^) Lagrangians are presented. The finite and gauge invariant sums of one-loop amplitudes are listed in the next 



three lines (M^±' = J2j M^f) 
sum of all amplitudes is given. 



The numbers 1, 2, 3 denote the row of diagrams on the Fig. H. In the last line the 



In the determination of D* D77 and B* S77 a sizable contribution from /?-like elec- 
tromagnetic terms (14) has been found |27|. Therefore we have to investigate their effect in the 
1)0 — > 77 decay amplitude. The considerations of Eq. ( [T^ ) gives us additional diagrams which are 
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^ ' — ® ® ^ ' ^=^(| 



Figure 4: The diagrams with one /3-hke (jlj) couphng (described by •), which give vanishing 
amphtudes. 



/ \ I \ 

II II ~K 

DO \ /' z?o /' ^0 

DO /' \ DO ;' \ DO ;' 

^ ' ® ® I J ® ® ^ ® g 

Figure 5: The diagrams which give nonzero amphtudes with one /3-hke couphng. 
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given in Fig. |^ and ^, where the /3 vertex is indicated by •. The nonzero parity violating parts of 
the one loop diagrams containing /? coupling are (Fig. |5|) 



.G 



1 



mz)— 4aie ga{j3 -\ 



1 



1 ml. 



./2 mc'(mz) + 2 A*) 167r2 3 



(-) 



2 ,(mz3)2 



/3.6 



V2 



G3(m., ^ + A^, H^) - G3(m., AJ, ^Hi)] }, 



(23) 



(24) 



while the parity conserving parts of the amplitudes arising from the one loop diagrams with /3 
coupling are 

Gf 



{+) 



/9.1 



=aie a(/3 + 



1 



1 



1 



V2m^ ' mc mo + 2A* 127r2 

aie2a(/3 + — ) 



1 



1 



/3.2 



M, 



{+) 



G 



rric {mo + 2A*) 127r 



(25) 



-M, 



(+) 



/3.1 



/3.3 



4=a,eV(/? + [/2(.«x, + A,) + /,(m^)] + 



+ 



VudVcdlhim^^mD + A^) + /i(m^)] |, 



(26) 



(27) 



{+) 



/3.4 



zaie ga{p + 



1 



1 



\/2m£) 



mc m£) + 2 A* Gvr^ 



■ 1 



Ks'K;! [2^1 ("^^) + ("^15 + A,)G3(mi^,mD + A„ -^)] + 
+ yudycd\j^h{m^) + {rriD + ^d)G^{m^, ruD + A^, -^)] }, 



(28) 



(+) 



G, 



/3.5 



^-^-,„^ + 2(A, - A*) 



1 



M, 



(+) 



/3.6 



y V* 

^US ^CS 



+Vudv:, 



Gf 2/^2 
— ==aie 5ra(/3 , 

(h^ + a:) 



mo + 2(Ad - A*; 



- [l2imK,mD + A,) - hiniK, ^ + A^)] + 



■ [l2{m^,mD + Ad) - /2(m^, ^ + A^)] }, 



(29) 



+ A, - A*)^^^""^' + 2 ; ' (rno + A, - A*) 



) + Jr^" ^ .. G,{mK,mn + A., 



)] 



+ 



™" ^ G3 (-., T + -T) + T^^^^^G, {m^,mn + A„ -^)] }. 



(!r^ + A,-A*) 



(^ + A, - A* 



(30) 



The amplitudes with /3 coupling are not finite and have to be regularized. We use the strict MS 
prescription A = 1 as in Q (note that in A = has been used). We take /i = 1 GeV ~ A-^ 

Hi. 



as m 
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In order to obtain the value of (5 we use the available experimental data from D*+ D'^j, 
£)*o _^ jjO^ decays. For instance, one can use the recently determined D*~^ decay width r{D*^) = 
96 ± 4 ± 22 keV || together with the branching ratio Br{D*+ D+-f) = (1.6 ± 0.4)% g. At 



tree level one has 



e2 



r{D*+ ^ D^j) = ^{IJ^^ - (31) 

2 

with = ^^^^(1 ^) the momentum of outgoing photon. Using the experimental data and 

rur = 1.4 GeV one arrives at P = 2.9 ± 0.4 GeV-i f\, where the errors reflect the experimental 
errors. 

On the other hand one can also use the ratio of partial decay width in D*^ system T{D*^ 
D^-f) : r(L>*° L>°7r°) = (38.1 ±2.9) : (61.9±2.9), where the experimental errors are considerably 
smaller than in the previous case. At tree level one has 

with k^,kn the momenta of outgoing photon and pion respectively. Using nic = 1-4 GeV, g = 0.59, 
f = = 132 MeV one arrives at (3 = 2.3 ± 0.2 GeV~\ [] where errors quoted again reflect 
experimental errors only. The /? coupling coming from from D*~^ ( |3l|) and D*^ ( |3^ ) are in fair 
agreement, but not equal. This signals that other contributions coming from chiral loops and higher 
order terms that would alter our determination of /5 might be important. Since the contribution of 
chiral loops to T(D*~^ D'^j) are approximately 50%, while for D*^ — > D'^j are about 20% |^6|| , 
we use in our numerical calculations the value of /5 = 2.3 GeV~^ obtained from T{D*^ D^l)- 
Results are shown in Table ^ using (3 = 2.3 GeV~^ and rric = 1.4 GeV. Inspection of Tables |l| and 
^ reveals that for the real parts of the amplitudes, the contributions of Figs. |2|, |^ and of Fig. ^ 
are comparable in size. However, the decay rate is dominated by the contribution of the imaginary 
part of the parity-violating amplitude, which arises from the one loop diagrams of Fig. ^. For the 
parity-conserving amplitude, the contributions of SD, anomaly and /3-like terms are comparable in 
magnitude. 

Due to the suppression of 02 in comparison to ai, we do not include diagrams proportional to 
a2 in the calculation of terms with 

Using short distance contributions, the finite one loop diagrams and the anomaly parts of the 
amplitudes (shown in Figs. ^ and with numerical values of the amplitudes as listed in Table 
one obtains 

5r(L»° ^ 77) = 1.0 X 10"^ (33) 

This result is slightly changed when one takes into account the terms dependent on (3 (0). The 
branching ratio obtained when we sum all contributions is 

Br{D^ -f-f) = 0.95 x 10"^. (34) 

By varying /? within 1 GeV~^ < /? < 5 GeV^^ and keeping g = 0.59 it 0.08, the branching ratio 
is changed by at most 10%. On the other hand, one has to keep in mind that the loop contributions 

*There is also a solution of (^ij) /3 — 0.09 ± 04 GeV~^ that, however, does not agree with the determination of {3 
from D*° decay. 

^The other solution is f3 = —3.6 ± 0.2 GeV~^ that does not agree with D*'^ data. 
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Diag. 


7 s 

M;^[x 10-10 GeV] 


/ 1 \ 

M^+)[x 10-10 GeV-i] 







-2.69 


/3.2 





2.69 


(3.3 





2.11 


PA 


0.88 


-0.007 


/3.5 





0.51 


/3.6 


-2.88 


-0.52 




-2.00 


2.09 



Table 2: Table of nonzero contributions of the amplitudes coming from the diagrams with /3 coupling (Fig. In 
the last line the sums of the contributions are presented. We use f} — 2.3 GeV"^, nic = 1.4 GeV. 



involving beta are not finite and have to be regulated. We have used MS scheme, with the divergent 
parts being absorbed by counterterms. The size of these is not known, so they might influence the 
error in our prediction of the branching ratio. Note also that changing a would affect the predicted 
branching ratio. For instance, if the chiral corrections do decrease the value of a by 30% this would 
decrease the predicted branching ratio down to 0.5 x IQ-^. 

4 Summary 

We have presented here a detailed calculation of the decay amplitude 77, which includes 

short distance and long distance contributions, by making use of the theoretical tool of Heavy 
Quark Chiral Perturbation Theory Lagrangian. Within this framework, the leading contributions 
are found to arise from the charged vr and K mesons running in the chiral loops, and are of order 
0{p'^). These terms are finite and contribute only to the parity violating part of the amplitude. 
The inclusion of terms of higher order in the chiral expansion is unfortunately plagued with the 
uncertainty caused by the lack of knowledge of the counterterms. As to the parity conserving part 
of the decay, it is given by terms coming from the short distance contribution, the anomaly and 
from loop terms containing the beta coupling, the latter giving most of the amplitude. The size of 
this part of the amplitude is approximately one order of magnitude smaller than the parity violating 
amplitude, thus contributing less than 20% to the decay rate. Therefore, our calculation predicts 
that the D — > 27 decay is mostly a parity violating transition. 

In addition to the uncertainties we have mentioned, there is the question of the suitability of 
the chiral expansion for the energy involved in this process; the size of the uncertainty related to 
this is difficult to estimate. Altogether, our estimate is that the total uncertainty is not larger than 
50%. Accordingly, we conclude that the predicted branching ratio is 

Br(Z)° ^ 77) = (1.0 ± 0.5) X 10-^ (35) 

The reasonability of this result can be deduced also from a comparison with the calculated decay 
rates for the — > p{uj)'y, which are found to be expected with a branching ratio of approximately 
10-6 |,0,|T1. 

We look forward to experimental attempts of detecting this decay. Our result suggests that the 
observation of D — > 27 at a rate which is an order of magnitude larger than (^) , could be a signal 
for the type of "new physics", which leads to sizable enhancement Q| of the short-distance c — > U7 
transition. 
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A List of chiral loop integrals 



Here we list dimensionally regularized integrals needed in evaluation of xPT and HQxPT one-loop 
graphs shown on Fig. & 



1 



1 



h{m) 



1 1 



with 



(27r)4-^ (g2_m2)(g.w- A) 167r2 A 



/2(m,A), 



/i(m) = m In I — ^ I — m A, 



/2(m, A) = -2A2ln 



m 



4A2f(^) +2A2(1 + A), 



(36) 
(37) 

(38) 
(39) 



where A = f — 7 + ln(47r) + 1 (in calculation A = 1), while F{x) is the function calculated by 
Stewart in |26|, valid for negative and positive values of the argument 



X 



X 
„2 



vr 

--tan 



-1 



In ( 3; + \/x2 — 1 



\x\ < 1 
Ixl > 1 . 



(40) 



The other integrals needed are (for k = 0) 

1 



1 1 



Go(m, A, v-k), 



(27r)4 (g2 _ m2)((5 + k^ - m^){q-v - A) Wir'^ vk 

Go{in, A, v-k) = h^{in, A) - h^{m, A + vk) - 7r[/i(m, A) - h{m, A + v-k)], 



(41) 



where 



/i(m, A) = < 



Arctan 
i In 



A 



-v/m2 — A2 
m 



\m\ > lAI 



A - VA2 - m? 



vr 



(42) 



and 



IfjL 



d^-'q 



+ sign(A)-; \m\ < |A| 



g^'''G^{m,A,vk) + 



(27r)4-^ (g2 _ m?){{g + kf - m'^){q-v - A) 167r2 

(w'^r + Fv^)G4(m, A, w A:) + k^'k'^G^im, A, A:) + ?;''v''G6(m, A, vk) 



(43) 
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with 



G3{m, A, v-k) 
G4{m, A, v-k) 
G^{m, A, v-k) 
GQ{m, A, v-k) 



m 

2v-k 
1 

v-k 

1 

v-k 
1 

2v-k 



Go{m, A, /c) 



1 



l2im,A) - l2{m,A + v-k) +A + ^, (44a) 



1 



4u-/c 

l2(m, A + V- /c) — G^im, A, f • /c) 



.2(A + 7;-A:) 

- ^iVo("i,0) + AG2(m, A,i;-A;) - Giim.A + v-k) 



him, A) - him, A + v-k) 



(44b) 

(44c) 
(44d) 



where him, A) is defined in (|37|) and NQim, k"^) in (^ 

In calculation we also need several other integrals that have been calculated for the 
mDV^", kl = kl = 0, v-ki = v-k2 = ^ 



d^q 1 _ 1 

(2^)4 iiq + A;i)2 _ m2)((g - yta)^ - m2)(^-t; - A) ~ Tfo? 
d^g 1 



Go(m, A,mz)), 
1 



(45) 



(27r)4 iiq + A:i)2 - m2)((g - yta)^ - m2)(g2 _ ^2)(^.„ _ a) 167r2 ^' ^^^^ 

where 

G'o(m, A,m£,) 



ArriD 



f - /i(m, A - i^^) Jm2 - (A - ^)2 - 



^-him,A + ^))Jm?-iA + ^)^-i5 



(47) 



2/i(m, 



i6 



Moim,A,mD) 



-^-^ I - 2h\m, ^) - 2/i2(m, A) + h\m, A - ^) + h\m, A + 



mo ■ 
2 . 



+ in In 



+ zvr In 



A - ^ - V"^' - (A - ^)" - i<5 

-A - ^ - - (A + T)' - ^'^ 

A + iVm'^ - A2 - i(5 



+ 



(48) 



-A + iV"i2 - A2 - i5 



with 5 > an infinitesimal positive parameter. 
The chiral loops needed are 



d^-'q 



(27r)4-^(g2_rn2)((g + A:)2-m2) IQn^ 
Noim, fc2) = -A + 1 - (^) + In K| , 



Noim,k^ 



(49) 



where 



Hia) = < 



2(1- v^4/a - 1 Arctan 



2 l-yn:4^-<^ In 



^1 - 4/a + 1 



yr^47^-i 



:0 < a < 4 



i7r0(a — 4) > ; otherwise 



(50) 
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and for kf = k2 = 



with 



q^'q'' 



1 



(27r)4-^ (g2 - m2)((g + kiY - m?){{q + fes)^ - m?) IGtt^ 

Mi[m, ki-k2 



g'"'M2im,krk2 



krk2 



M2{m,kvk2) 



ki-k2 



M^{m, ki-k2 



Mi{m,krk2) = 

4 la 



+ 1 - ^~Arcth 



Arcth 



+ '-(u2(^- 

1 



+ Li2 



Li2 



1 + ^ 



+ Li2 



(51) 



(52a) 

(52b) 
(52c) 



where we have abbreviated a = 2ki-k2/'m? and 
arithmic function. 



y^l + 2m? /ki-k2, while Li2(x) is a polylog- 
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